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We have found the complete spectrum and eigenstates for 
harmonic oscillations of ideal spherical and cylindrical shells, 
both being infinetely thin. The spectrum of the cylindrical 
shell has an infinite number of Goldstone modes correspond- 
ing to folding deformations. This infrared catastrophe is over- 
come by accounting for curvature-dependent part of energy. 
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INTRODUCTION 

The oscillating surfaces have numerous applications in many areas of 
physics, spanning from biophysics and chemical physics to high energy physics. 
1,2 ' 3 In this article, we study the oscillation spectra of spherical and cylindri- 
cal surfaces, considering them as infinitely thin elastic membranes. We also 
calculate their heat capacities. Corresponding calculated values of heat capac- 
ity are valid at sufficiently low temperature T « <d D , where Qd is Debye 
temperature, otherwise the effects of discreteness become essential. To our 
knowledge and big surprise, this problem has not yet been solved, despite of a 
vast literature on the subject. 4 



THE SPECTRUM AND EIGENSTATES OF SPHERICAL AND 
CYLINDRICAL SHELLS 
A. Spherical Membrane 

The Lagrangian of a spherical membrane reads : 



where R is radius of spherical membrane, dQ = sin 9d6d(j) and, // and A are the 
Lame coefficients. u r , ue, are components of the displacement vector and 
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Uee, U 4,4,1 U00 are components of the strain tensor. The latter are determined 

by relative change of distance between two points at a deformation : ds" 2 = 
ds 2 {l + U a pdx a dxf3) and can be expressed in terms of the displacement vector 
and its derivatives as follows: 5 

TJ - 1 dUe + Ur (?) 

Uee ~ R~de + r (2) 

tt 1 du <t> , Ud +a , Ur /q\ 

u ** = R^eW + R cotd+ R (3) 

1 ,du* nx 1 dua 

2^ = 5 (^-«.*co t ») + — - (4) 

Equations of motions for oscillations with the frequency u read : 
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- [X + 3^ + 2(A + /i)/Cj + + {pu?R 2 + 2 M ) We = (6) 

'\ + fi + 2(A + p)/c) (^gg) + (A + 3^ + 2(A + „)/c) (gj^) 

+ (A + 2^ + 2(A + „)/c) (-^ ) + + 

+ (pw 2 /? 2 + 2/i)«^ = (7) 



where /C = • There exists a symmetric solution (the "breathing" 

mode), uq = u§ = , u r = const ^ 0. The frequency of this solution is u s = 



4 ^2^ ■ From other solutions, we first consider those independent on <p. There 
exist three branches of such oscillations. The first two (" longitudinal" modes) 
correspond to = 0, Uq ^ 0, u r ^ 0. The frequency of these oscillations are 
labeled by an integer I : 

[4A + 2fi + (A + 2fi)l(l + 1)] 
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V /(4A + 2^ + (A + 2fi)l(l + 1))2 - 16^(A + /x)(/(Z + 1) - 2) 

± w 2 {) 

In the third branch (a "shear" mode) uq = , u r = and ^ 0. Its 
frequencies are 

^ = -^(/(/ + l)-2) (9) 

where / = 2, 3, ... Explicit expressions for the displacements for the longitudinal 
branches are: 6 

dP t ( cos 6) . . 

U(f> = u e = C l \ ' (10) 

ctu 

where C is constant amplitudes and P;(cos#) are the Legendre polynomials. 
The displacements in the shear branch solution read : 



u r = ue = = C 



dPi(cos 6) 
d9 
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In order to get eigenstates for m ^ 0, we employ an infinitesimal rotation 
around x axis which turns an arbitrary vector r into : 



r = 0(uj)r = 1 + tog x r 



(13) 



Here — ig x is the generator of the rotation around x axis. 7 It generates a 
transformation of the displacement field u{r). 

u \r) — 0(a;)'u(0~ 1 (c<j)r) = u{r) — ug x u(f) + u(g x f).Vu(r)... (14) 

If u (f) was a solution of the system (5-7), then u (r) is also a solution. The 
operation we performed to find new solutions is readily recognizable as the Lee 
translation introduced in theory of fiber bundles. 8 
After a simple transformation, we find an operator Q, which, acting on a 



solution fixed by two components 
the same frequency. 



( \ 



\ u <j> ) 



, generates another solution with 



Q = 



sin (j)-^ + cot 9 sin 0^ 



sin 6 



V 



cos 4> 
sin 6 



sin J| + cot 9 sin <f>^r 



(15) 
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The operator Q is the analog of the operator of total angular momentum J in 
quantum mechanics, but it is reduced to a subspace of tangential components 
of the displacement vector u. The Lagrangian (1) and equations of motion 
(6-7) are invariant with respect to the transformation (15). Acting by the 
operator Q onto solutions with m = (10-12), we generate solutions with 
all possible values of m . They belong to the same eigen-frequencies as the 
initial solutions with m — 0. The displacements for solutions belonging to the 
longitudinal branches with m/0 read : 



u { r ] (x) = e im<t> i m 



mx(l - x z )^P, (m) (x) - (1 - x z )^P, 



X 
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u f\ x ) = _ e im ^ m+1 m(l - x 2 )^Pr\x) 



(m) , 



(17) 



The displacements for oscillations of the shear branch with m^O are : 



u^\x) = - e m H m+1 m(l - x 2 )^Pr\x) 



(m), 
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u f\ x ) = e irn fi r - 



mx{l - x 2 )^P ; (m) (x) - (1 - x 2 )^P l 



X 



(19) 



Given the known frequency spectrum ujj, the thermal energy can be cal- 



culated as U 



° ~~ \dT'V 



hijOj 

exp^r- — l 



and the heat capacity can be calculated as 



As always, low-lying rotational frequencies u>^ ot ^ = are smaller than 

low-lying oscillatory frequencies by a factor ~ \f%% where a is characteristic 
distances between atoms ( discreetness effect ) and m is the mass of electron. 
That means that at temperature T » Tr = -^jp, the rotational degrees of 
freedom become substantially classical until k < y^T. Their contribution to 
the heat capacity is Cr = ^£(3). Elastic oscillations remain frozen till a 
larger temperature T ~ >> T R . 

B. Cylindrical shells 

In order to find the spectrum and corresponding eigenstates of a cylindrical 
membrane, we follow the same steps as we did for spherical shell. For cylin- 
drical case, the components of strain tensor in terms of cylindrical coordinates 
r , 0, z read: 5 

tj 1 du 4> , u r T j du z 1 du z d U4> 

The Lagrangian reads : 

' du r \ 2 / du^ \ 2 / du z \ 2 



- ± J Rd^dz{U zz + U H ) 2 -tij Rd^dz{U 2 zz + 2U% + U 2 ^) (21) 
where R is the radius of the cylinder. Equations of motions are : 



u r = (A + 2fi)IC^ + iCRX^ (22) 
o<p oz 

u t = -(A + 2,)^ - (A + „)*|* + -f^((A + 2,)^ - 

(23) 

where K = pR2a)2 ^ A+2 ^- ) . Specifying the z dependence of all displacement 
components as a plane wave e tkz , and their angular dependence as e m ^, we 
find a system of linear equations for the amplitudes u r , ue , u^. Its secular 
equation reads : 

{u 2 - to 2 ) (u 6 - (uj 2 {R 2 k 2 + m 2 ) + u 2 (R 2 k 2 + m 2 + l))u; 4 

'{{R 2 k 2 + m 2 ) 2 + 5R 2 k 2 + m 2 ) - ^\R 2 k 2 ^ 2 

+4 Wt Vt 2 -" a 2 )#* 4 ) =0 (25) 

where u; 2 = and a; 2 = One of its solutions (symmetric) is u 2 = uj 2 s . 
For k — 0, the remaining cubic equation can be solved explicitly for arbitrary 
m. The three solutions are : 

lo 2 = 0, ui 2 = (m 2 + IV 2 , lo 2 = m 2 uj 2 (26) 



+ [u*u 2 J 



The first of them implies that there exists infinite number of Goldstone modes 
corresponding to the folding deformations. For m — 0, the cubic equation 
gives the following solutions : 



The physical reason of the peculiarity in the spectrum of cylindrical shell 
is an opportunity to fold a cylinder conserving its straight-lines and its cross- 
section contour length. Such a deformation is strainless. The infinite number 
of folding modes leads to divergence of the heat capacity at any temperature. 
This difficulty can be avoided by introducing terms with higher derivatives into 
the Lagrangian. Geometrically, it means that the Gaussian curvature and the 
mean curvature must be included. In more physical or chemical terms it means 
that not only the variation of distances between atoms, but also variation of 
valence angles must be accounted for. Below we write down explicit formu- 
lae for the Gaussian curvature K and the mean curvature M for cylindrical 
membranes in linear approximation in displacements. 



2 r 2 d2 2 

co = k R u t 



(27) 




l(tfk 2 + i) ± yW^ 2 + 1)) 2 - 



uj 2 s - uj 2 )R 2 k 2 



(28) 
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K = 



1 d 2 u r 



(29) 
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1 u r 1 9V <9V ,„n\ 



The term in Lagrangian depending on curvature reads : 



AL 



-I 



^(K-K ) 2 + d -(M-M ) 2 



d 2 S (31) 



where K = and M = ^ are the values of K and M for undeformed 
cylinder. With this term in Lagrangian we find the gaps A m in the spectrum 
corresponding to the former folding modes: 



A m = y^(m 2 -1) (32) 

Note that the gaps in spectrum due to curvature effect are propotional to R~ 2 
and, therefore, are small in comparison to gaps due to strain effects (oc Z? -1 ). 
By perturbation method it is possible to find small k - dependent corrections 
to these frequencies : Su 2 = [2k 2 R 2 (m 2 — 1)] — ^ or u m(k) = 
A 2 m + V 2 k 2 , where V 2 = ^{2R 2 (m 2 - 1)] - ( "^gf )2 - A clear tendency 
is that the band becomes more flat with the growth of m. 
For cylindrical shell, we can write the thermal energy in the following form : 

U=J^f ^# (33) 

In the temperature interval %A 2 « T « huj t , we can neglect all " elas- 
tic branches" and substitute the summation over the "valence" branches by 
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integration. It leads to the energy 

_ °- 601 

and the heat capacity : 

_ 1.803 



\ PR 2 



-T 3 (34) 



(35) 



h 2 >| ^(o* - 2«?) 
For T « hA 2 , we can neglect the sum over m, since only A 2 and V2 contribute 
to the energy. The internal energy becomes exponentially small : 
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U = V2jrhT-^e-— (36) 

and the heat capacity : 

CONCLUSIONS 

We have found complete spectrum of oscillations for spherical and cylindrical 
membranes and corresponding eigenstates. We have showed that cylindrical 
shell has an infinite number of Goldstone modes due to folding deformations 
which do not change the elastic energy of shell. We have also showed that 
this catastrophe can be avoided by introducing curvature-dependent terms 
to the Hamiltonian. In addition, heat capacities of spherical and cylindrical 
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shells are calculated. They can be measured experimentally as an additional 
contribution to the heat capacity of a solvent in a diluted solution of long 
molecules. 
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